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A new global regression analysis method is discussed that predicts wind 
tunnel model weight corrections for strain gage balance loads during a wind 
tunnel test. The method determines corrections by combining “wind on” model 
attitude measurements with least squares estimates of the model weight and 
center of gravity coordinates that are obtained from “wind off” data points. 

The method treats the least squares fit of the model weight separate from the 
fit of the center of gravity coordinates. Therefore, it performs two fits of “wind 
off” data points and uses the least squares estimator of the model weight as 
an input for the fit of the center of gravity coordinates. Explicit equations 
for the least squares estimators of the weight and center of gravity coordinates 
are derived that simplify the implementation of the method in the data system 
software of a wind tunnel. In addition, recommendations for sets of “wind off” 
data points are made that take typical model support system constraints into 
account. Explicit equations of the confidence intervals on the model weight and 
center of gravity coordinates and two different error analyses of the model weight 
prediction are also discussed in the appendices of the paper. 
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= matrix containing regressors of the first least squares problem 

= matrix containing regressors of the second least squares problem 

= axial force of the alternate balance axis system (see Fig. 5 or Ref. [1], Fig. 1) 

= weight correction to the axial force of the alternate balance axis system (see Fig. 5) 
= axial force that is caused by aerodynamic loads on the wind tunnel model 
= axial force that acts in the primary balance axis system (see Fig. 2) 

= axial force that is caused by the weight of the wind tunnel model 
= axial force that is caused by the weight of the sleeve (metric part of the balance) 

= axial force that is caused by the combined weight of the model and the sleeve 
= vector containing the dependent variables of the first least squares problem 
= vector containing the dependent variables of the second least squares problem 
= vector that describes the location of the center of gravity of the wind tunnel model 
= components of vector D in the primary balance axis system 
= least squares estimators of d x ,d y , d z 

= set of unit basis vectors of the primary or alternate balance axis system 
= gravity vector in the primary balance axis system 

= components of the gravity vector G in the primary balance axis system 
= gravitational acceleration (equals 9.81 [m/sec 2 ] or 32.17 [ft/sec 2 ] at sea level) 

= index of an observation 
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= number of regressors used in Ref. [4] 

= auxiliary matrix used by the first least squares problem 
= auxiliary matrix used by the second least squares problem 
= auxiliary vector used by the first least squares problem 
= auxiliary vector used by the second least squares problem 
= combined mass of the wind tunnel model and the metric part of the balance 
= residual mean square 
= total number of observations 

= normal force of the alternate balance axis system (see Fig. 5 or Ref. [1], Fig. 1) 

= weight correction to the normal force of the alternate balance axis system (see Fig. 5) 

= normal force that is caused by aerodynamic loads on the wind tunnel model 
= normal force that acts in the primary balance axis system (see Fig. 2) 

= normal force that is caused by the weight of the wind tunnel model 
= normal force that is caused by the weight of the sleeve (metric part of the balance) 

= normal force that is caused by the combined weight of the model and the sleeve 
= total number of regression model parameter 

= pitching moment of the alternate balance axis system (see Fig. 5 or Ref. [1], Fig. 1) 

= weight correction to the pitching moment of the alternate balance axis system (see Fig. 5) 
= pitching moment that is caused by aerodynamic loads on the wind tunnel model 
= pitching moment that acts in the primary balance axis system (Fig. 2) 

= pitching moment that is caused by the weight of the wind tunnel model 
= pitching moment that is caused by the weight of the sleeve (metric part of the balance) 

= pitching moment that is caused by the combined weight of the model and the sleeve 
= total number of data points used to fit the weight and center of gravity coordinates 
= rolling moment of the alternate balance axis system (see Fig. 5 or Ref. [1], Fig. 1) 

= weight correction to the rolling moment of the alternate balance axis system (see Fig. 5) 

= rolling moment that is caused by aerodynamic loads on the wind tunnel model 
= rolling moment that acts in the primary balance axis system (Fig. 2) 

= rolling moment that is caused by the weight of the wind tunnel model 
= rolling moment that is caused by the weight of the sleeve (metric part of the balance) 

= rolling moment that is caused by the combined weight of the model and the sleeve 
= side force of the alternate balance axis system (see Fig. 5 or Ref. [1], Fig. 1) 

= weight correction to the side force of the alternate balance axis system (see Fig. 5) 

= side force that is caused by aerodynamic loads on the wind tunnel model 
= side force that acts in the primary balance axis system (Fig. 2) 

= side force that is caused by the weight of the wind tunnel model 
= side force that is caused by the weight of the sleeve (metric part of the balance) 

= side force that is caused by the combined weight of the model and the sleeve 
= percentage point of the f-distribution (see, e.g., Ref. [4], p. 515, Table A. 3) 

= weight vector, i.e., the force of gravity on a body 

= length of the weight vector W ; equals the “scalar” weight of the wind tunnel model 
= least squares estimator of w 
= weight estimate of a single data point 

= coordinate that defines the direction of the axial force -or- regressor variable 
= generic matrix that holds the regressor values of a multivariate least squares problem 
= vector containing the coefficients of the first least squares problem 
= vector containing the coefficients of the second least squares problem 
= coordinate that defines the direction of the side force -or- response variable 
= vector that holds the response values of a multivariate least squares problem 
= yawing moment of the alternate balance axis system (see Fig. 5 or Ref. [1], Fig. 1) 

= weight correction to the yawing moment of the alternate balance axis system (see Fig. 5) 
= yawing moment that is caused by aerodynamic loads on the wind tunnel model 
= yawing moment that acts in the primary balance axis system (Fig. 2) 
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= yawing moment that is caused by the weight of the wind tunnel model 

= yawing moment that is caused by the weight of the sleeve (metric part of the balance) 

= yawing moment that is caused by the combined weight of the model and the sleeve 

= coordinate that defines the direction of the normal force 
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= parameter that defines the percentage points of the t-distribution (see Ref. [4], p. 515) 
= generic vector that holds the regression coefficients of a least squares problem 
= generic regression coefficients 
= least squares estimators of /3 i,/?2,/33 

= dimensionless unit gravity vector in the primary balance axis system (Fig. 2) 

= components of the unit gravity vector T in the primary balance axis system (Fig. 2) 

= prediction error estimate of the axial force of the balance 
= prediction error estimate of the normal force of the balance 
= prediction error estimate of the model weight 

= angle between two vectors in the primary balance axis system (Fig. 2) 

= pitch angle of the wind tunnel model 
= index of a data point 
= standard deviation 


I. Introduction 


A strain-gage balance is used during a wind tunnel test to measure forces and moments that act on a 
wind tunnel model. These loads can be determined, for example, in three steps. First, the raw gage outputs 
of the balance are measured as absolute voltages when the loads act on the wind tunnel model. Then, the 
natural zeros, i.e. , the known gage outputs of the absolute load datum of the balance, are subtracted from 
the raw outputs to get the output changes relative to the absolute load datum of the balance. Finally, the 
output changes are used as input for a load iteration scheme that computes the balance loads (see Ref. [1] 
and [2] for a detailed description of the load iteration scheme). 

The user of an “unpowered” wind tunnel model is primarily interested in aerodynamic loads that act on 
the model in a “wind-on” test condition. Unfortunately, a balance cannot directly distinguish between the 
part of its electrical outputs that is caused by aerodynamic loads and the part that is caused by the weight 
of the model (see Fig. 1). Therefore, the predicted balance load, i.e., the result of the load iterations, is 
always the sum of an aerodynamic load and a load that is caused by the weight of the model. This situation 
can be expressed by using the following relationships for the forces and moments of a six-component balance 
(the subscript “ bal ” marks the total predicted load; the subscript “aero” marks the part that is caused by 
aerodynamic loads; the subscript “ weight ” marks the part that is caused by the total model weight): 
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It was mentioned earlier that the change of the strain-gage outputs relative to the natural zeros of the 
balance is used as input for the iterative balance load calculation because the natural zeros represent the 
outputs of the balance gages in an assumed “weightless” condition. Therefore, balance loads associated with 
weight effects are not just caused by the weight of the wind tunnel model. The weight of the balance sleeve, 
i.e., of the metric part of the balance, also has to be considered for the sake of completeness even though 
its contribution is small when compared with the weight of a typical wind tunnel model. Consequently, the 
total weight effects on the balance load measurements of a “wind-on” data point can be expressed by the 
following equations (the subscript “ weight ” marks the total model weight; the subscript “ model ” marks the 
part of the total model weight that is exclusively caused by the wind tunnel model; the subscript “sZeeue” 
marks the part of the total model weight that is caused by the metric part of the balance): 
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The term “model weight” is frequently used in subsequent sections of the paper. Its usage implicitly 
assumes that the weight of the metric part of the balance is a part of the model weight. 

Equations (la) to (If) indicate that loads caused by the weight of the model, i.e., A we i g ht , ■ ■ • , Y we i g ht , 
have to be subtracted from the computed “raw” balance loads, i.e., At, a i, TbaZ, so that the desired 
aerodynamic loads, i.e., A aero , Y aero , are obtained. Consequently, a process has to be defined that 
determines balance loads caused by the total model weight for every “wind-on” data point. The development 
of such a process starts with a “first principles” analysis of all loads that are associated with the gravitational 
acceleration (i.e., A we i g ht , ■■■, Y we i g ht for a six-component balance). This analysis and the subsequent 
definition of loads associated with the model weight are described in the next section of the paper. 


II. Definition of Model Weight Corrections 


It is necessary to describe all balance load components that are associated with the gravitational ac- 
celeration so that wind tunnel model weight corrections can be predicted. First, a balance axis system 
needs to be defined before the loads can be described. It is assumed that (i) the balance axis system is a 
right-handed coordinate system and that (ii) both force and moment vectors follow the sign convention of 
the related Cartesian coordinates. Figure 2 shows the corresponding coordinate system. It is also called the 
“primary” balance axis system. This axis system and its load sign definitions are the basis for all derivations 
that are discussed in the paper. It also has the useful property that the cross-product operator of vector 
algebra can rigorously be applied during the derivation of the model weight corrections. At this point of the 
analysis it is important to explicitly mention that the signs of two of the six loads used in the “primary” 
balance axis system do not agree with load sign definitions that are commonly adopted in experimental aero- 
nautics. This difference of the sign definitions, however, can easily be overcome by performing two simple 
load transformations (for more details see the discussion given in Section V of the present paper). 

The balance loads associated with the gravitational acceleration are determined by (i) the combined 
mass of the wind tunnel model and the metric part of the balance and (ii) the location of the gravity vector 
in the balance axis system. The gravity vector G can be described using the components of the gravitational 
acceleration in the balance axis system. Then, we get: 
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The length of the gravity vector equals the gravitational acceleration itself at the given elevation of an 
observer. Therefore, we know that 


|G| = g = yj gl + gl + g 2 z 

where, for example, the gravitational acceleration at sea level equals 

g = 9.81 {m/sec 2 } ...or .. . g = 32.17 [ft/ sec 2 ] 


(36) 


(3c) 


It is convenient to describe the gravity vector as a dimensionless quantity because the current analysis 
is primarily interested in a description of the direction of the gravity vector in the balance axis system. 
Consequently, the following equation is valid: 
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The length of the dimensionless gravity vector T is one because the following relationship is valid: 
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In the next step it is necessary to connect the total mass of the wind tunnel model and the dimensionless 
unit gravity vector T to the weight vector W (see Ref. [3], p. 127, for a definition of the weight vector). We 
know from Newtonian Mechanics that 


W = m ■ G = m • g ■ V 
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Then, after introducing the model weight as a scalar quantity, we can write: 

w = m • g 

Consequently, the weight vector can be expressed as follows: 


lx 


1 

3 

l 

ly 

= 

ly ■ w 

lz. 


lz ■ W_ 


(7) 


(8) 


(9«) 


The components of the weight vector are identical with the axial force, side force, and normal force that 
are caused by the weight of the model. Therefore, we get: 
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Finally, after equating the right-hand side of Eq. (9 6) with the right-hand side of Eq. (9a), we get for the 
forces that are caused by the model weight in the “primary” balance axis system the following relationships: 


FORCES CAUSED BY THE MODEL WEIGHT 

A weight 

= 7 x-w 

(10 a) 

S weight 

= 7 y‘W 

(106) 

weight 

= 7 z-W 

(10c) 


The moments associated with the weight of the model can be obtained using a similar approach. These 
moments have to be described relative to the balance moment center as the center of gravity of the model 
is not necessarily identical with the balance moment center. The present analysis simply assumes that the 
origin of the balance axis system is the balance moment center. Then, the center of gravity of the model can 
be described by the using the following vector (see also Fig. 3): 


D = 



( 11 ) 


Vector D describes the model’s lever arms in the balance axis system. The torque caused by the model’s 
weight with respect to the balance moment center is defined by the cross-product of the “lever arm” vector 
D and the “weight” vector W (see also Ref. [3], p. 302). In addition, we know that the components of the 
torque are the rolling, pitching, and yawing moment. Therefore, we can write: 
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The cross-product defined in Eq. (126) needs to be computed. In general, the cross-product of two 
vectors can be described as a determinate. Then, we get: 
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The forces A we i g ht , S we i g ht, and N we i g ht used in Eq. (12c) can be replaced by the right-hand sides of 
Eqs. (10a), (106), and (10c). Then, Eq. (12c) becomes: 
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In the last step, after equating the right-hand side of Eq. (13) with the right-hand side of Eq. (12a) 
we get for the moments in the “primary” balance axis system the following relationships: 
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MOMENTS CAUSED BY THE MODEL WEIGHT 
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The six equations (10a), (106), (10c), (14a), (146), and (14c) define all forces and moments that are 
caused by the model weight. These loads are known for any “wind-on” data point if the following three sets 
of parameters are given: Set 1 => the components ixilyilz of the unit gravity vector of the “wind-on” 
data point in the balance axis system; Set 2 =>■ the weight w of the wind tunnel model; Set 3 => the 
center of gravity coordinates d x ,d y ,d z of the wind tunnel model in the balance axis system. 

The components of the unit gravity vector change for each “wind-on” data point because any model 
attitude change also changes the components of the unit gravity vector in the balance axis system. Therefore, 
lx-, lyi Iz have to be measured for each “wind-on” data point in “real-time” by using the wind tunnel’s 
model attitude measurement process. The weight of the wind tunnel model ( w ) and its center of gravity 
coordinates ( d x , d y , d z ), on the other hand, are four physical constants of the chosen model configuration as 
long as (i) no hardware is added to or removed from the model and (ii) the model’s weight distribution is not 
changed. These four constants have to be determined experimentally for a “real-world” wind tunnel model 
because the complexity of its geometry and weight distribution makes an analytic determination impractical. 

The four constants may be obtained, for example, by applying a regression analysis, i.e., a least squares 
fit, to a series of “wind-off” data points. This approach has two major advantages. First, only a regression 
analysis in combination with a suitable experimental design and a tight control over the most important 
measurement error sources can ensure that both systematic and random errors contained in numerical 
estimates of the model’s weight and center of gravity coordinates are within a wind tunnel customer’s 
tolerances. In addition, aerodynamic loads of “wind-off” data points are zero by design (see also Fig. 4). 
Consequently, the measured strain-gage outputs and predicted loads of “wind-off” data points are exclusively 
caused by model weight effects. This observation can be expressed by using the following set of equations 
(the subscript “ bal ” marks the total load; the subscript “aero” marks the part that is caused by aerodynamic 
loads; the subscript “ weight ” marks the part that is caused by the total model weight): 


“WIND OFF” TEST CONDITION 
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The calculation of least squares estimators of the weight and the center of gravity coordinates of the 
wind tunnel model is described in more detail in the next section of the paper. 
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III. Global Regression Analysis 


A. General Remarks 

A new method is derived in this section that may be used to determine least squares estimators of the 
weight (w) and center of gravity coordinates (d x , d y , d z ) of a wind tunnel model by exclusively using explicit 
equations. This novel approach greatly simplifies the implementation of the method in the data system of a 
wind tunnel facility as complexities associated with the coding of a global regression analysis algorithm are 
avoided. In principle, the method performs a global regression analysis of a set of “wind-off” data points. It 
makes the assumption that the least squares fit of the weight of the model can be separated from the least 
squares fit of the center of gravity coordinates. In other words, the new method performs two “coupled” 
least squares fits of “wind-off” data points as it uses the result of the first least squares fit, i.e. , the least 
squares estimator of the model weight, as an input for the second least squares fit. 

The global regression analysis of the set of “wind-off” data points is explained in detail in the next 
sections of the paper. More general information about the theoretical background of the global regression 
analysis approach can be found in the literature (see, e.g., Refs. [4], [5], [6], and [7] for more details about 
global regression analysis). 


B. Least Squares Estimator of the Model Weight 

First, the new method fits the model weight by exclusively using the three force components that are 
measured for each “wind-off” data point. Index £ describes an individual “wind-off” data point assuming 
that a total number of q data points were recorded. Then, after applying Eqs. (10a) to (10c), we get: 


; 1 < £ < Q 

1 < £ < q 

; i < e < q 

Equations (16a) to (16c) can be summarized in vector format. Then, we get: 
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Now, a more simple notation of the above two vectors is introduced. In addition, the model weight w 
is replaced by its least squares estimator w. Then, we get: 
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Then, Eq. (17) can be expressed as follows for a single data point with index £: 
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Consequently, we get the following equation if Eq. (19a) is expanded to all “wind-off” data points: 
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For convenience, the two vectors associated with all data points can be abbreviated as follows: 
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Then, after using Eqs. (20a) and (206) to simplify Eq. (196), the equations of the first global regression 
analysis problem of the new method are obtained. We get: 


First Global Regression Analysis Problem ==^ A 39X1 • X lxl = B 39X1 (21) 

The solution of the first global regression problem, i.e., the least squares estimator of the model weight, 
is simply given by the following matrix equation: 


LEAST SQUARES ESTIMATOR OF THE MODEL WEIGHT 






/ — T — \ _1 


X lxl 
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( A ' A ) - 

X 7 lxl 
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The right-hand side of Eq. (22) is the new method’s least squares estimator of the model weight. It 
is possible to derive an explicit formula for the model weight after recognizing that the right-hand side of 
Eq. (22) is simply the product of two scalars. Therefore, by inspection, we can write: 

[ 7x(0 + 7 y(0 + 7*(0 ] = 

v ✓ 

see also Eq. (6) 

[ 7x(0 • A bal ( 0 + 7 „(0 • S bal (0 + 7,(0 • N bal (0 } (236) 

Then, after using the right- hand sides of Eq. (23a) and (236) in Eq. (22) and after recognizing that the 
“inverse” of the “matrix” used on the left- hand side of Eq. (23a) equals 1/q, we get: 

w = • Y [ 7x(0 • Aai(0 + 7 yiO'ShaliQ + 7z (0 ' N hal (0 ] (24) 

q «=1 

It is possible to simplify Eq. (24) even further. The summation term on the right-hand side of Eq. (24) 
can be written as a scalar product of two vectors. Then, we get: 

r(0-w(0 = 7*(0-^M0 + ■ s b ai(0 + (25) 


1.0 = q (23 a) 

«=i 



9 

American Institute of Aeronautics and Astronautics 


In general, the scalar product of two vectors can also be expressed as the product of their lengths and 
the cosine of the angle between them. Therefore, we can write: 

m ■ W(0 = I r(0 I • I W(0 I • COS T 7 (0 (26) 

The length of the two vectors used in Eq. (26) is given by the following equations: 

| m I = 1-0 (27a) 

|W(0| = yj + S ^) + NUO (27 b) 

In addition, we conclude from Eq. (7) that the two vectors are parallel to each other. Therefore, we get: 

Eq. (7) =► T(0 II W(0 =► »J(0 = 0 =► cos 77(0 = 1 (27c) 

Now, after using the right-hand sides of Eqs. (27 a), (27 b) and (27c) on the right-hand side of Eq. (26), 
we get for the scalar product of the two vectors the expression: 

r (0 • w (0 = J Aitf) + s 2 bal ( o + nu o (28) 

Then, after comparing Eq. (28) with Eq. (25), we get the following identity: 

7x(0 -4,01(0 + 71/(0- SbaliO + 7*(0- = yj + S 6o«(0 + NUO (29) 

Finally, after substituting the summation term of Eq. (24) by the right-hand side of Eq. (29), we get the 
following explicit solution for the least squares estimator of the model weight (the corresponding confidence 
interval on the model weight is derived in App. 1): 


EXPLICIT SOLUTION OF THE MODEL WEIGHT 

™ = jr ■ E ] / A L® + SL(0 + NU 0 (30) 

q «=i 


A useful interpretation of Eq. (30) can be made => the least squares estimator of the model weight 
equals the arithmetic mean of the individual weights that are calculated for each “wind-off” data point. 
It is also concluded from Eq. (30) that the accuracy of the model weight prediction is directly linked to 
the accuracy of the axial, side, and normal force predictions of the balance. The accuracies of these force 
predictions, assuming they are described in engineering units, differ significantly for a typical six-component 
strain-gage balance because the axial force capacity is often one order of magnitude smaller than the normal 
force capacity. In other words - the prediction error of the axial force is expected to be substantially smaller 
than the corresponding error of the normal force. This conclusion suggests that the wind tunnel model 
should be pitched up as high as possible in order to maximize the contribution of the axial force in Eq. (30). 
Unfortunately, the pitch angle range of a typical wind tunnel model goes from p. « —5° to p ~ +20°. 
Consequently, the error associated with the normal force prediction will dominate the total error associated 
with the model weight prediction for the majority of wind tunnel models (an upper bound of the model 
weight prediction error is derived in App. 2; a more general error analysis is described in App. 5). 

C. Least Squares Estimators of the Center of Gravity Coordinates 

Now, the remaining least squares estimators of the center of gravity coordinates of the wind tunnel 
model have to be determined. The new method uses a second least squares fit of the “wind-off” data points 
for this purpose. The second least squares problem is defined by assuming that the weight of the model is 
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known as the result of the first least squares fit. This weight, i.e. , the least squares estimator w, is defined in 
Eq. (30) . In addition, the second least squares fit is constructed by exclusively using the predicted moments 
of the “wind-off” data points. Then, after applying Eqs. (14a) to (14c), we get the following equations for 
a single “wind-off” data point with index £ : 


Rbal (0 

= 7 z(0 

• w 

unknown 

7v(0 

• w 

unknown 

Pbal(ti) 

= 7*(£) 

• w 

- 

unknown 

7 *(£) ■ 

■ w • 

unknown 

Ybal(0 

= 7 y{0 ■ 

w • 

unknown 

7*(£) ' 

w • 

unknown 


Equations (31a) to (31c) can be expressed in matrix format. Then, we get: 


1 < £ < q 

1 < £ < q 
i < £ < q 


- 7 z(0 ■ w 

+ 7 y{Q-W 


+ 7z(0 • w 

0 

- 7®(0 • w 


- 7 y(0 ■ W 
+ 7z (£) • W 
0 



dx 


Rbal(€) 


dy 

= 

Pbal(0 


d z 


_Y b aM)_ 


model attitude measurement and weight from Eq. (30) 


unknown 


balance load 


(31a) 

(316) 

(31c) 


(32) 


In the next step, abbreviations of the matrix and vectors used in Eq. (32) are introduced. In addition, 
the coordinates d x , d y , and d z are replaced by their least squares estimators d x , d y , and d z . Then, we get: 


K(£) 


3X3 


0 + 7z(0 ■ w 

-7 z($)-w 0 

+ 7 y(0 ■ W - 7x(0 • W 


L(0a xl 


X 


3X1 


Rbal (£ ) 
Pbal(0 

Y ba l(0 

d x 

dy 

dz 


- 7 y(0-w 
+ 7x(0 • w 

0 


(33a) 


(336) 


(33c) 


Consequently, Eq. (32) can be expressed as follows for a single “wind-off” data point with index £: 


K(<£) 


• X 

3X3 ^3X1 


= L(0a xl 


(34a) 


In addition, we get the following equation if Eq. (34a) is expanded to all “wind-off” data points: 


'K(1) 3x3 ‘ 


'LWaxa' 

K(2) 3x3 

• x 3xl = 

£(2) 3x1 

- ^(^3X3 - 


- ^( 9 ) 3 x 1 - 


(346) 


For convenience, the matrix and vector associated with the data points can be abbreviated as follows: 


A 


3 q X 3 


'K(1) 3x3 ' 

K(2) 3x3 

-K(g) 3x3 _ 


(35a) 
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(356) 


B, 


L(l)sx 1 

£( 2) 3x1 


Mq) axi. 

Then, after using Eqs. (35a) and (356) to simplify Eq. (346), the equations of the second global regression 
analysis problem are obtained. We get: 


Second Global Regression Analysis Problem 


X, 


= B, 


(36) 


The solution of the second global regression problem can be described by the following matrix equations: 


LEAST SQUARES ESTIMATORS OF THE CENTER OF GRAVITY COORDINATES 





( AT 

A ) 

7 3X3 

dx 

• X 3xl = (a t b) 

(37a) 

3X1 





Y 

7V 3X1 

dy 

d z 

= V A ) - 

x 7 3X3 

A T B (376) 

' 7 3X1 


The solution of the second global regression problem is expressed in two different ways. The first matrix 
equation, i.e., Eq. (37 a), represents the “normal equations” of the global regression problem. They define 
the solution as a 3 x 3 linear system of equations that has to be solved for the three unknown center of 
gravity coordinates contained in vector X. The second matrix equation, i.e., Eq. (376), is an “academic” 
solution of the linear system that is defined in Eq. (37a). It expresses the solution as a matrix multiplication 
assuming that the inverse matrix used on the right-hand side of Eq. (376) has been determined. 

An explicit solution of Eq. (37a) is derived in App. 3. In addition, the confidence intervals on the center of 
gravity coordinates are derived in App. 4. It is always helpful if the center of gravity coordinate predictions 
can independently be verified. The center of gravity of a model with symmetric weight distribution, for 
example, has to be located in the plane that is defined by the axial force and the normal force in the balance 
axis system. Therefore, the new method’s predicted least squares estimator d y of the center of gravity 
coordinate d y has to be close to zero for a model with a symmetric weight distribution. 

IV. Selection of “Wind— Off” Data Points 

The “wind-off” data points for the global regression analysis of the weight and center of gravity coor- 
dinates of the wind tunnel model have to be selected carefully so that (i) the confidence intervals on the 
corresponding least squares estimators are meaningful and that (ii) the linear system associated with the 
second global regression analysis problem, i.e., Eq. (37a), has a unique solution. It is shown in the literature 
that a unique solution of a general global regression analysis problem exists if the columns of the matrix of 
the regressors, i.e., matrix A of Eq. (35a), are linearly independent (for more detail see Ref. [4], p. 69). 

Different approaches are used in the aerospace testing community to select “wind-off” data points 
for the determination of the model’s weight and center of gravity coordinates. The authors, for example, 
recommend to use at least five “evenly” distributed data points for the global regression analysis in order 
to have sufficiently small confidence intervals on the least squares estimators (see App. 1 and App. 4 for a 
discussion of the connection between the size of the confidence interval on the least squares estimators and the 
number of “wind-off” data points). The bounds of the data point set could be defined by the model’s largest 
positive pitch angle (fj, m ax) and the model’s largest negative pitch angle (fimin)- In addition, the “wind-off” 
data point at “zero” pitch angle (/i = 0°) could be used as a “confirmation point” in order to independently 
validate the prediction of the model weight and center of gravity coordinates that are the result of the least 
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squares fits of the five “wind-off” data points. Table 1 below lists this first set of recommended “wind-off' 
data points including the confirmation point: 


Table 1: First set of “wind-off” data points ("f y = 0 ; Sbai = 0 ; q = 5 ; one confirmation point). 



Pitch Angle 

Roll Angle 

-A-bal 

Nbal 

1 

A^(l) = ^ 0 

0° 

M ,az(l) 

AM 1) 

2 

M(2) — f^min 0.25 • (/imax ft min ) 

0° 

A-bal (2) 

AM2) 

3 

M(3) = ftmin H" 0.50 * (ftmax fumin') 

0° 

Abal(3) 

Nbal (3) 

4 

M(4) — ftmin 0.75 • (//max fumin') 

0° 

Abal (4) 

Nbal (4) 

5 

= ftmax 0 

0° 

Abal (h) 

AM5) 

6 

/i(6) = 0° (Conf. Pt.) 

0° 

\Abal\ ~ 0 

I Nbal I ~ w 


It was mentioned during the discussion of the explicit solution of the model weight, i.e. , Eq. (30), that 
the accuracy of the normal force prediction has a substantial influence on the total error associated with 
the model weight calculation (see also the error analysis that is discussed in App. 2). Therefore, whenever 
possible, the authors recommend to record “wind-off” data points with the model in both “upright” and 
“inverted” orientation, i.e., for the model roll angles 0° and 180°. In this case, the normal force will have 
both positive and negative signs which may help reduce a potential “hidden” bias in the experimental data 
that may exist if the normal force of the “wind-off” data points never changes sign. Table 2 below shows, 
for example, “wind-off” data points for a model in both “upright” and “inverted” orientation (the two 
“wind-off” confirmation points at “zero” pitch angle, i.e., fi = 0°, are also listed). 


Table 2: Second set of “wind-off” data points (q y = 0 ; Sbai = 0 ; q = 10 ; two confirmation points). 


£ 

Pitch Angle 

Roll Angle 

A-bal 

Nbal 

1 

Mi) 

0° 

Abal (1) 

N ba l (1) 

2 

M2) 

0° 

Abal( 2) 

Nbal(2) 

3 

M3) 

0° 

Abal (3) 

Nbal(3) 

4 

M4) 

0° 

Abal( 4) 

NbaM) 

5 

M5) 

0° 

Abal (5) 

AM5) 

6 

M 6) = M 1 ) 

180° 

Abal( 6) ~ —A b al( 1) 

Nbal (6) » -Nbal(l) 

7 

M 7) = M2) 

180° 

Abal (7) ~ —Abal (2) 

N b al (?) W ~Nbal(2) 

8 

M8) = M3) 

180° 

Abal (8) ~ —Abal (3) 

Nbal (8) » ~Nbal(3) 

9 

M9) = M 4 ) 

180° 

Abal (9) ~ —Abal (4) 

Nbal (9) « ~Nbal(4) 

10 

Mio) = m(5) 

180° 

Abal (10) ~ — Abal (5) 

Nbal(10) « -N ba l(5) 

11 

Ml) = 0° (Conf. Pt. 1) 

0° 

\Abal\ ~ 0 

| N ba l | ~ W 

12 

M 2 ) = 0° (Conf. Pt. 2) 

180° 

\Abal \ ~ 0 

A^feaJ S3 W 


The selection of “wind-off” data point suggested in both Table 1 and Table 2 may potentially be 
improved by choosing data points such that the variance in the regression coefficients, i.e, the variance in 
the model weight and in the center of gravity coordinates, is minimized. This goal can be achieved by 
using a point selection that maximizes the determinates of A A and A T A (a discussion of the related 
experimental design problem can be found in Ref. [4], p. 505/506). 

Unfortunately, the aerospace testing community traditionally uses two load sign definitions that do 
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not follow the signs of the related Cartesian coordinates of a right-handed coordinate system (see, e.g., 
the balance axis system definition used in Ref. [1], p. 3). Therefore, load transformations may have to be 
performed in order to apply the new method to traditional wind tunnel test data. These transformations 
are discussed in the next section of the paper. 


V. Balance Load Transformations 


The “primary” balance axis system is shown in Fig. 2. It is a right-handed coordinate system. It 
also uses sign definitions for both the forces and moments that match sign definitions of related Cartesian 
coordinates. Only in this case the cross-product operator of vector algebra can be used for the calculation 
of moments that are associated with model weight effects. Therefore, the “primary” balance axis system 
was chosen for the rigorous derivations of all model weight corrections that are discussed in this paper. 

The aerospace testing community, on the other hand, uses an “alternate” balance axis system. This 
“alternate” balance axis system is depicted in Fig. 5. Its coordinate and load sign definitions were directly 
copied from the literature (taken from Ref. [1], p. 3, Fig. 1). Unfortunately, the “alternate” balance axis sys- 
tem has two load components that contradict the sign definitions of related Cartesian coordinates. Therefore, 
the “alternate” balance axis system has the disadvantage that the cross-product operator of vector algebra 
cannot be used for the calculation of moments that are associated with model weight effects. Consequently, 
the “alternate” balance axis system is not suited for a rigorous “hrst-principles-based” derivation of the 
model weight corrections. 

Differences between the “primary” and “alternate” balance axis systems and associated load sign def- 
initions can most easily be recognized by viewing Fig. 2 and Fig. 5 side-by-side. It can be seen that the 
coordinate axes and load sign definitions of Fig. 2 and Fig. 5 are identical with the exception of the sign 
definitions of the axial and normal forces. The positive axial force of the “primary” balance axis system 
points in the direction of the positive x-axis. The positive axial force of the “alternate” balance axis system 
points in the direction of the negative z-axis. The positive normal force of the “primary” balance axis 
system points in the direction of the positive z-axis. The positive normal force of the “alternate” balance 
axis system points in the direction of the negative z-axis. 

Fortunately, it is possible to use all results presented in this paper for the “primary” axis system with 
the traditional load sign definition of the aerospace testing community. It is only required to perform two 
load transformations. First, assuming that balance loads follow the traditional sign convention used in the 
aerospace testing community (see Fig. 5), the axial force and normal force of the “wind-off” data points are 
transformed from the “alternate” to the “primary” balance axis system by using the following relationships: 


■A-bal ~ 

= (-1) • AF 

(38 a) 

Nbal = 

-- (-1) • NF 

(38 b) 

All other balance load components defined in the “alternate” 
with the results presented in this paper. Therefore, we get: 

balance axis system can directly be used 

^bal 

= SF 

(38c) 

Fbal 

= RM 

(38 d) 

Fbal 

= PM 

(38e) 

Ybal 

= YM 

(38/) 


Then, the model weight corrections of the six balance loads are determined as outlined in the previous 
sections. Finally, the model weight corrections of the axial and normal force have to be transformed back 
from the “primary” to the “alternate” balance axis system. The following transformations apply: 


AF weight 

= (- 1 ) ' 

Ayjeight 

(39a) 

N F We igbit 

= (- 1 ) 

N weight 

(39 b) 
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All other model weight corrections, i.e., S weigh t , Rweight , Pweight , and Y weight can directly be applied 
to the original balance loads as no sign change is required. The following relationships apply: 


R R weight 

^weight 

(39 c) 

RM weight 

Rweight 

(39 d) 

PMweight 

Rweight 

(39e) 

Y Myjeight 

Y we igJit 

(39/) 


It is important to point out that the Cartesian coordinate definitions of the “primary” and “alternate” 
balance axis systems are identical. Therefore, no sign transformations need to be applied to the least squares 
estimators of the center of gravity coordinates. They are valid in both balance axis systems. 

VI. Conclusions 

A new method for the prediction of wind tunnel model weight corrections for strain-gage balance loads 
was introduced. The new method became possible after rigorously applying the ideas of global regression 
analysis to a set of “wind-off” data points that are needed for the prediction of the least squares estimators 
of the weight and center of gravity coordinates of the wind tunnel model. Figure 6 summarizes basic elements 
of the new method to predict model weight corrections. 

For the first time, a complete set of explicit equations is published that may be used for the determination 
of the least squares estimators of the model weight and center of gravity coordinates. In addition, explicit 
equations of the confidence intervals on the model weight and center of gravity coordinates are listed in the 
appendices of the paper. This complete set of explicit equations greatly simplifies the use of the new method 
in the data system of a wind tunnel facility as complexities associated with the implementation of a global 
regression analysis algorithm are avoided. 

The new method uses two separate least squares fits for the determination of the model weight and 
center of gravity coordinates. First, it is proven that the least squares estimator of the model weight equals 
the arithmetic mean of the model weights of the individual “wind-off” data points that are used for the 
global regression analysis. Then, least squares estimators of the center of gravity coordinates of the model 
are derived as the explicit solution of a 3 x 3 linear system of equations. In addition, two different error 
analyses of the model weight are discussed in the appendices of the paper. Finally, recommendations for 
suitable sets of “wind-off” data points are made that may be used for the determination of wind tunnel 
model weight corrections. 
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Appendix 1: Confidence Interval on the Model Weight 


The wind tunnel model weight was introduced as a regression coefficient in an earlier section of the 
paper. The model weight was obtained after fitting strain-gage balance forces as a function of related 
model attitude measurements. The reporting of results of such a regression analysis becomes more complete 
if confidence intervals on the regression coefficients are estimated. Therefore, explicit equations for the 
confidence interval on the model weight are discussed in this section. These equations were derived from 
more general relationships that are reported in the literature (see Ref. [4], pp. 41-42). This approach was 
selected as a derivation of the general equations of the confidence intervals goes beyond the scope of the 
present paper. First, the general equations of the confidence interval on a regression coefficient of a “no- 
intercept” model are developed using information that is provided in the literature. Then, it is shown how 
these general equations can be applied to the wind tunnel model weight. 

The discussion of the general equations of the confidence interval on a regression coefficient starts with 
a “no-intercept” regression model as this type of model was used to obtain the model weight (see Eqs. (17), 
(196)). This “no-intercept” model can be described as follows (see Ref. [4], Eq. (2.48)): 

V = Pi ■ x (40) 

where x is the regressor variable, y is the response variable, and pi is the regression coefficient. Furthermore, 
it is assumed that a total of n observations, i.e. , (yi,Xi), i = 1,2, ...,n, are to be used to determine the 
least squares estimator Pi of the slope Pi. Then, Eq. (40) can be expressed for each observation by using 
the following equation: 

yi = pi ■ Xi ; 1 < i < n (41) 

Similarities between the general “no-intercept” model discussed in Ref. [4], pp. 41-42, and the model 
used for the determination of the model weight exist that can be recognized after comparing terms of Eq. (196) 
with terms of Eq. (41). Table 3 below summarizes these similarities: 

Table 3: Comparison of the Nomenclature used for the Model Weight Determination. 


Reference [4] 

present paper 

Comments 

Pi 

W 

regression coefficient; model weight 

Pi 

W 

least squares estimator; Eq. (30) 

n 

3 q 

total number of “observations” 

•&1 1 %2 1 • • • J ^TL 

7a(l)> 7y(l), 7z(l)> ••■.7 7 z(q) 

regressor variable values 

Vl, 2/2, • • • , Vn 

Abalip ), Sbal( 1 ), N ba l( 1 ), • • • , S ba l (?) , A T ba l(q) 

response variable values 


The analytic solution of the regression coefficient of the “no-intercept” regression model is given in the 
literature as follows (from Ref. [4], Eq. (2.50)): 


Pi 


^ Vi Xi 


y x i 

i= 1 


(42) 


Now, the question comes up if the regression solution of the model weight (i.e., Eq. (24) and Eq. (30)) 
can be obtained “by inspection” from the general solution of the regression coefficient of the “no-intercept” 
model that is given in Ref. [4] (see also Eq. (42) above). In that case, the confidence interval on the model 
weight could directly be deduced from the confidence interval on the regression coefficient of the general 
“no-intercept” model that is listed in Ref. [4]. The answer to the question can be found by using the 
nomenclature comparison shown in Table 3 for the replacement of pi, Xi, and yi in Eq. (42). Then, we get: 
- = Aal{ l)'7x(l) + S bal (l) ■ 7^(1) + ••• + s ba i(q) ■ 7y(q) + N ba i(q) ■ r y z (q) 

W 7x(l) + 7y (1) + ••• + l 2 y{q) + 7z (?) 1 j 
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Equation (43) can also be expressed as follows: 


E [ A bal{£.) ~ 7x(£) + SbaliQ' 7j/(0 + Nbal{£) ■ 7z(£) ] 
«=1 

W = 9 

£ [ 7*(0 + 7jJ(0 + 7z(0 ] 

€=i 

We also know from Eq. (23a) that the following relationship applies: 




q 

E 

£=i 


7x(0 


7, 2 (0 


7?(6 


(44) 


(45) 


Then, after replacing the denominator of Eq. (44) with the left hand side of Eq. (45) and after rearranging 
the factors of the products of the summation terms of the numerator of Eq. (44), we get: 



(46) 


As expected, the right hand side of Eq. (46) exactly matches the right hand side of Eq. (24). Therefore, 
the general least squares solution of the regression coefficient of the “no-intercept” model from the literature 
(i.e., Eq. (42) above, or, Ref. [4], Eq. (2.50)) matches the solution for the model weight that was developed 
in an earlier section of the current paper (i.e., Eq. (24) and Eq. (30)). Consequently, it is possible to use the 
formula of the confidence interval from the literature (Ref. [4], Eq. (2.53)) to directly develop the confidence 
interval on the model weight. It is only necessary to replace the generic nomenclature used for the definition 
of the confidence interval in Ref. [4] with the nomenclature that is associated with the determination of the 
model weight. 

Reference [4] constructs the confidence interval on the regression coefficient Pi by making the normality 
assumption on the errors of the least squares fit. Then, the 100(1 — a) percent confidence interval on the 
slope /3i of the “no-intercept” model is defined as follows (from Ref. [4], Eqs. (2.52), (2.53)): 


CONFIDENCE INTERVAL ON THE SLOPE OF A “NO INTERCEPT” MODEL 
(. . . taken from Ref. [4], Eqs. (2.52), (2.53), . . .) 


Pi - t c 


/ 2 , 


\ 


MSr„ 


< Pi < Pi + t 


E 

2=1 


/ 2, n- 1 


\ 


MS Res 


E x i 


(47a) 


2=1 


where the residual mean square MS Res is defined as 


MS Res 



n n 

E y i ~ & E yi Xi 

. 2=1 2=1 


(476) 


Equations (47a) and (476) can be used to construct the confidence interval on the model weight. It 
is only required to replace the generic terms by the corresponding wind tunnel model specific terms. For 
example, the following two relationships can be obtained from Table 3 and Eq. (46) above: 

n = 3q (48) 
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w 


1 

9 


(49) 


01 


■ E [ 7® (0 • Abalit) + 7 y(0 ' S bal(0 + 7z (0 ' N bal (£) ] 
5=1 


It was shown above that the right hand side of Eq. (24) equals the right hand side of Eq. (46) (i.e., 
Eq. (49)). In addition, it was proven earlier that the right hand side of Eq. (24) equals the right hand side 
of Eq. (30). Therefore, after replacing the right hand side of Eq. (49) with the right hand side of Eq. (30), 
we get the following simplified relationship: 

01 = = V • i V + NUZ) (50) 

q 5=1 

Now, after using the contents of Table 3 to replace Xi and after applying Eq. (23a) to the result, we get: 

n q 

E x i = 7*(1) + 7y(l) + ••• + 7* ( 9 ) = X![7x(£) + % l ^ + 7^(0] = 9 (51) 

2=1 £=1 

Similarly, after using the contents of Table 3 to replace we get: 

n q 

E y i = ^LU) + S 6al(!) + ••• + N lM) = E[ A tal(0 + SU0 + N U0] (52) 

2 = 1 f=l 

We also know, after rearranging Eq. (42), that the following relationship applies: 

n n 

E y* Xi = 0i ■ x i ( 53 ) 

i = 1 i = 1 

The right hand side of Eq. (53) can be expressed differently. First, the right hand side of Eq. (50) can 
be used to replace /3i . Then, the right hand side of Eq. (51) can be used to replace the summation term on 
the right hand side of Eq. (53). Consequently, we get the following relationship: 


n 

E y* Xi 

2=1 


E \j A bal(Q + S bal(0 + N bal(£) 
5=1 


(54) 


At this point, all generic terms used in Eqs. (47a) and (476) (i.e., terms associated with /3±, n, Xi , and yt) 
can be replaced with corresponding wind tunnel model specific terms (i.e., terms associated with w , q, Ab a i, 
Sbai, an d Nbai) by using the right hand sides of Eqs. (48), (50), (51), (52), and (54). Then, after simplifying 
and rearranging the result, we get the following final explicit equations for the confidence interval on the 
model weight: 
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EXPLICIT SOLUTION OF THE CONFIDENCE 
INTERVAL ON THE MODEL WEIGHT 

W - t a/2 3 g — , • C(q) < w < w + t a/2 3 g_ x • C{q) (55 a) 


where the weight w 


is given by Eq. (30) and the parameter C\q) is defined as 


C{q) 


E{q) - F{q) 
Q ■ [ 3 <? - 1 ] 


(556) 


and the auxiliary parameters E(q) and F(q) are defined as 


E(q) 

F(q) 


Y, [ ^L;(0 + SbaliO + ^bal(0 ] 

S=1 

r i 2 

' y y/AUt) + S bal(0 + N bal(£) 
q €=1 


(55 c) 


(55d) 


One parameter in Eq. (55a) still needs to be specified before the confidence interval on the model weight 
can explicitly be calculated. It is the percentage point of the f-distribution it a/2 3? _ 1 ). This constant 
depends on (i) the confidence interval choice (e.g., 95 % or 99 % confidence interval) and (ii) the number of 
“wind-off” data points (i.e., q) that are used for the least squares fit associated with the model weight. The 
limiting case, i.e., the use of a single data point for the least squares fit, must be avoided. This case will 
not allow for the calculation of the confidence interval on the model weight as the confidence interval will be 
zero. This assertion can easily be proven “by inspection” because the following relationships apply: 

q = 1 => E{ 1) = F(l) => C{ 1) = 0 (56) 

The confidence interval is zero because the factor C(l) is zero. Therefore, the number of data points 
needs to be greater than one in order to obtain a meaningful non-zero value for the confidence interval. 
Table 4 below lists the percentage points for different confidence interval choices and data point numbers. 

Table 4: Percentage Points of the t Distribution for p = 1 (from Ref. [4], p. 515, Tabic A. 3). 


Cl 

a 

a/2 

q 

II 

GO 

kQ 

1 

l— 1 

^ a/2 , 3q — l 

95 % 

0.05 

0.025 

3 

8 

2.306 

95 % 

0.05 

0.025 

5 

14 

2.145 

95 % 

0.05 

0.025 

10 

29 

2.045 

99 % 

0.01 

0.005 

3 

8 

3.355 

99% 

0.01 

0.005 

5 

14 

2.977 

99% 

0.01 

0.005 

10 

29 

2.756 


Table 4 shows that already a total number of five data points (i.e., q = 5) leads to percentage points of 
the f-distribution that are near the value of 2 for the 95 % confidence interval and near the value of 3 for 
the 99 % confidence interval. 
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Appendix 2: Upper Bound of the Model Weight Prediction Error 

It was shown in Eq. (30) that the least squares estimator of the model weight equals the arithmetic 
mean of the weights that are computed for each “wind-off” data point. Therefore, it is possible to perform 
an error analysis of the model weight prediction that (i) assesses contributions of the different balance load 
measurements to the total weight prediction error and (ii) defines an upper bound of the prediction error. 
The analysis of the model weight prediction error starts by expressing Eq. (30) in a slightly different format. 
The least squares estimator of the model weight can be described by the following equation: 

™ • X) where W '(£) = \! A bal(0 + S baM) + N baM) ( 57 ) 

q t=l 

The model’s roll angle is often either 0° or 180° when “wind-off” data points are recorded. Then, the 
gravity vector’s component in the side force direction of the balance axis system is zero. Therefore, we get: 

l v « 0 => S bal « 0 (58) 

Now, after applying Eq. (58) in Eq. (57) and after dropping the “wind-off” data point index £ for 
simplicity, the summation term u/(£) defined in Eq. (57) becomes: 

= yj A Li + *L (59) 

Equation (59) assumes that both the axial and normal force measurements are “perfect.” In reality, 
however, the prediction of balance loads is always associated with errors. These errors are primarily caused by 
load magnitude measurement and load alignment errors that occurred during the calibration of the balance. 
Therefore, the summation term defined in Eq. (59) has to be replaced by the sum of the “true” model weight 
and an associated weight prediction error. We get: 

w’ => w’ + Aw' (60) 

Similarly, the balance loads need to be replaced by the sum of the unknown “true” load and the 
associated load prediction error. Then, we get for the axial and normal forces the relationships: 

Abai => A ba i + A A bai (61a) 

N ba i => N ba i + A N ba i (61 6) 

Consequently, after replacing the loads in Eq. (59) by the right hand sides of Eqs. (61a) and (615), the 
sum of the true model weight and the weight prediction error can be written as follows: 


w' + Aw' — \J [ A ba i + A A ba i } + [ N ba i + AN ba i } (62) 

Then, after expanding the brackets on the right hand side of Eq. (62), we get: 

w' + Aw' = \J [ Al al + N£ al } + [ A Al al + A Ng al ] + 2 • [ A bai ■ A A ba i + N ba i ■ A N ba i } (63) 

We also know, after comparing the order of magnitude of terms on the right hand side of Eq. (63), that 
the following inequality and approximation are valid: 


2- 

[ A ba i ■ AAbal , + N ba i 

■ AAbaZ ] 


» 


+ AA bal y 

(64a) 


1st order 

1st order 



2nd order 

2nd order 



second order terms 

= =► AA 

2 

bal 

+ 

AN Li * 

0 

(646) 
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In addition, the following abbreviation can be introduced: 


first order terms =4> e = 2 • [ A ba i ■ A A bal + N ba i ■ AN ba i ] 

Then, after using the relationships given in Eqs. (64a), (646), and (64c), Eq. (63) becomes 

/ml 


w 


Aw' 


, t— N 2 

ly bal 


(64c) 


(65) 


It is possible to simplify the right hand side of Eq. (65) by using a Taylor Series. In general, the Taylor 
Series approximation of a function in the neighborhood of an independent variable ip is defined as follows 


/(/ + 


/(/) 


fW 

1! 


/"(/) 

2 ! 


( 66 ) 


where e is a small value. Then, after comparing the right hand side of Eq. (65) with the left hand side of 
Eq. (66), we know that 

/(/ + e) = yj [ A% al + Ng al ] + e (67) 

Therefore, we need to find a Taylor Series expansion for the following function type: 


/(/; ) = y/lp~ 

where the first derivative of /(/) is given by the following expression: 


/'MO = 


1 

\P¥ 


(68a) 


(686) 


Then, after inserting Eqs. (68a) and (686) into the right hand side of Eq. (66) and after neglecting higher 
order terms by truncating the Taylor Series after the first derivative, we get: 


\[W 


\P¥ 


\P¥ 


We also know that 


/ — Alai + N? 


bal 


(69) 


(70) 


Then, after replacing if in Eq. (69) by the right hand side of Eq. (70) and after comparing the result 
with the left hand side of Eq. (65), we get: 


Aw' 


/AI 


bal ""H ^ bal 


{k. 


bal _r N ba i 


1 


V Al a i + N bal 


We also know, after rearranging Eq. (64c), that the following relationship is valid: 

-j- = A ba l ■ A A ba l + N b al ■ A N bai 

Then, after using the right hand side of Eq. (72) to replace the factor e/2 in Eq. (71), we get: 

A ba i ■ A A bai + N ba i ■ A N ba i 


Aw' 


A b al + N bal 


equals w' 


\fA, 


bal Al bal 


(71) 


(72) 


(73) 


equals w' 


We also know that the two square roots on the right hand side of Eq. (73) can be replaced by the left 
hand side of Eq. (59). Then, after simplifying the result, we get: 


w 


Aw' 


/ . 
w + 


Abal ■ A Abal + Nbal ‘ A N ba l 


w' 


(74) 
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Now, a term can be identified in Eq. (74) that may be used to study the model weight prediction error. 
We get the following estimate of this error after simply comparing the left and right hand sides of Eq. (74): 

Abai ■ A A ba i + N ba i ■ A N ba i 


A w' 


w' 


(75) 


It is possible to define an upper bound of the weight prediction error as the following inequality applies: 

A ba i ■ AA ba i + N ba i ■ AN bai ^ \A ba i\ ■ AA ba i + \N ba i\ ■ AN ba i 


Aw' 


< 


(76) 


w' w' 

In addition, the following trigonometric relationships between the axial force, normal force, weight 
estimate, and the pitch angle of the model are valid: 

\A ba i\ « w' • | sin n | (77 a) 

|A bai | « w' • | cos /x | (77 b) 

Finally, after using the right hand sides of Eqs. (77a) and (77 b) to replace the absolute values of the 
forces on the right hand side of Eq. (76), we get the following inequality for the upper bound for the error: 


UPPER BOUND OF THE MODEL WEIGHT PREDICTION ERROR 


Aw' 


AA, 


bal 


sin \i 


AN, 


bal 


COS [1 


(78) 


axial force contribution 


normal force contribution 


In general, the normal force capacity of most 6-component strain-gage balances is one order of mag- 
nitude larger than the axial force capacity. Then, assuming that errors are given in engineering units, the 
error associated with the normal force prediction is expected to be roughly one order of magnitude larger 
than the error associated with the axial force prediction because the electrical output range of a strain-gage 
is typically sized to match the related load range. This conclusion can be summarized as follows: 

A N bal « 10 • A A bal => A N ba i » AA ba i (79) 

The pitch angle of the model also has an influence on the weight prediction error as trigonometric 
functions of the pitch angle are used in inequality (78). Table 5 below compares, for example, the sine and 
cosine values for the typical pitch angle range of a wind tunnel model. 

Table 5: Sine and cosine values for a typical pitch angle range. 



sin /i 

COS (1 

0° 

0.000 

1.000 

5° 

0.087 

0.996 

10° 

0.174 

0.985 

20° 

0.342 

0.940 


Therefore, after combining inequality (79) with the expected values of the sine and cosine functions for 
the typical pitch angle range of a wind tunnel model, we get the following inequality: 

A N bai • | cos fj, | > AA ba i ■ | sin /x | (80) 

Finally, after comparing inequality (80) with inequality (78), it is concluded that the error associated 
with the normal force prediction, i.e. , A N ba i, is by far the largest contributor to the error associated with 
the model weight prediction over the typical pitch angle range of a wind tunnel model. This important 
conclusion can be illustrated by a numerical example. Let us assume that (i) the model is pitched to 20° 
and that (ii) the normal force error AN ba i is ten times as large as the axial force error AA ba [. Then, after 
evaluating all terms on the right hand side of inequality (78) and after simplifying the result, we see that 
the normal force contribution is 96.5 % of the upper bound of the model weight prediction error. 
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Appendix 3: Explicit Solution of the Center of Gravity Coordinates 


An explicit solution of the center of gravity coordinates can be obtained after applying Cramer’s Rule 
to the 3x3 linear system of equations that is defined in Eq. (37 a) (see, e.g., Ref. [5], pp. 323-324, for a 
description of Cramer’s Rule). First, the coefficients of matrix A T A and the components of vector A T B 
have to be determined by inspection. The following coefficients of matrix A T A are obtained: 


A t A = 


Mil U \2 Mi 3 
M21 M22 M23 

«31 u 32 W33 


where «n, M 12 , • • •, M33 are defined as 

Mn = 

U22 = 

U33 = 

U12 = 

U13 = 

U23 — 


follows: 


9 


w 2 ■ 

E 

Q 

[ 7y(C) 2 

+ 7,(C) 2 ] 

w 2 ■ 

E 

£=1 

Q 

[ 7 ,(0 2 

+ 7x(£) 2 ] 

w 2 • 

E 

«=i 


+ 7 y(Z) 2 ] 

9 

M21 


-w 2 ■ 

E [7x(0 '7,(0 
£=1 

9 

M31 


-w 2 • 

E [7,(0 '7,(0 

5=1 

Q 

M32 

= 

-w 2 • 

E [7,(0 ‘7,(0 

5=1 


Similarly, the following components are obtained for vector A T B: 


(81) 

(82a) 
(82 b) 
(82 c) 
(82 d) 
(82e) 
(82/) 


A 1 B 


^1 

v 2 

V3 


(83) 


where Mi, z> 2 , and M 3 are defined as follows: 


9 


Ml 

= w 

•E 

£=i 

[ Ybal{£) ' 

' 7,(0 

- Pbal (£ ) ' 

' 7,(0 ] 

(84a) 

v 2 

= VJ 

<? 

■E 

[Rbalti) 

• 7 ,(0 

— YbaliO 

• 7,(0 ] 

(84 b) 

V3 

= w • 

Q 

E 

[ Pbal(£,) ' 

7,(0 

~ Rbaltt) 

• 7,(0 ] 

(84c) 


f=i 


The coefficients Mu, M 12 , . . ., M 33 and the components Mi, M 2 , M 3 describe the left and right hand sides 
of the 3x3 linear system of equations that is defined in Eq. (37a). All these values have the factor w in 
common. Therefore, they can be divided by w without changing the solution of the linear system, i.e., of 
the center of gravity coordinates. Then, the coefficients Mu, . . ., M33 are replaced by rn, . . ., ^ 3 : 

rn = w ■ [7y(6 2 + 7z( 6 2 ] (85a) 

5=i 
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r 22 


w ■ 


(85 6) 


[ 7 ,(o 2 + 700 2 ] 

5=1 

^33 = W ■ [ 7x( 6 2 + 7,(0 2 ] (85c) 

5=1 

9 

ri2 = 721 = - w ■ [ 7x(0 ’7,(0 ] (85d) 

5=i 

9 

ri3 = r 3 1 = - w ■ 1 7,(0 ■ 7x(0 ] (85e) 

5=1 

9 

r 23 = C 32 = - w ■ ^2 [7,(0 ’7,(0] (85/) 

5=1 

Similarly, the components v 3l v 2 , v 3 are replaced by p 3 , p 2 , p 3 : 


9 


Pi 

ii 

[ Ybal(0 

■ 7 ,(0 

- Pbal (0 ' 

' 7,(0 ] 

(86a) 

P2 

ii 

[ RbaliO 

• 7,(0 

- Y bal (0 ■ 

• 7*(0 1 

(86 6) 

P3 

9 Wi 

H 

[ Pbal(0 ■ 

7x(0 

- Rbal(0 

• 7,(0 ] 

(86c) 


Finally, after applying Cramer’s Rule , we get an explicit solution of the center of gravity coordinates: 


EXPLICIT SOLUTION 

OF THE CENTER OF GRAVITY COORDINATES 


d x — 

Ri / 

A 

(87 a) 


dy = 

R2 / 

I? 

(87 6 ) 


d z = 

Ra / 

A 

(87c) 

where the four 

determinates R. 

Aii 

A 2 , and A 3 are defined as 



cn 

r 12 

c 13 


R 

= 

C 21 

722 

723 

(88a) 



1*31 

732 

C 33 




Pi 

n 2 

n 3 


Ri 

= 

P2 

1*22 

723 

(886) 



P3 

732 

733 




I’ll 

Pi 

C 13 


R'i 

= 

721 

P2 

723 

(88c) 



r-Ai 

Pa 

733 




I’ll 

n 2 

Pi 


Ra 

= 

721 

722 

P2 

(88 d) 



ra i 

732 

P3 
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Appendix 4: Confidence Intervals on the Center of Gravity Coordinates 


The center of gravity coordinates of the wind tunnel model were introduced as regression coefficients 
in an earlier section of the paper. These coefficients were obtained after fitting the predicted moments of 
a strain-gage balance as a function of (i) related model attitude measurements and (ii) the least squares 
estimator of the model weight. The reporting of the center of gravity coordinates becomes more complete 
if confidence intervals on the coordinates are estimated. Therefore, explicit equations for the confidence 
intervals on the center of gravity coordinates are discussed in this section. These equations were derived 
from more general relationships that are reported in the literature (see Ref. [4], pp. 93-94). This approach 
was selected as a derivation of the general equations of the confidence intervals goes beyond the scope of the 
present paper. First, the general equations of the confidence interval on a regression coefficient are developed 
using information that is provided in the literature. Then, it is shown how these general equations of the 
confidence intervals can be applied to the center of gravity coordinates. For simplicity, it is also assumed 
that the least squares estimator of the model weight, i.e. , w, can be treated like a constant in the regression 
model of the center of gravity coordinates. 

The discussion of the general equations of the confidence interval on a regression coefficient starts with 
a generic “multiple” linear regression model that is given in the literature. This model can be described as 
follows (from Ref. [4], Eq. (3.2)): 

y = p 0 + /3i • xi + p 2 ■ x 2 + ■■■ + Pk ■ x k + e (89a) 

where x is the regressor variable, y is the response variable, /Jo, • • ■ , Pk are the regression coefficients, and e 
is a small quantity. We know from Eq. (32) that the regression model associated with the center of gravity 
coordinates has no intercept (i.e., f3 0 = 0), and consists of three regression coefficients (i.e., k = 3). Therefore, 
after also dropping the small quantity e, Eq. (89a) can be written as follows: 

V = Pi ■ x-! + P 2 • x 2 + Pz • x 3 (89 6) 

Furthermore, it is assumed that a total of n observations, i.e., (jji, Xi), i = 1,2, ... ,n, are to be used to 

determine the least squares estimators Pi, p 2 , Pz of the coefficients Pi, p 2 , Pz of Eq. (896). Then, Eq. (896) 
can be expressed for each observation with index i by using the following equation (see also Ref. [4], Eq. (3.8)): 

Vi = Pi ■ xn + p 2 ■ x i2 + Pz ■ x i3 ; 1 < i < n (90) 

Reference [4] constructs the confidence interval on the three regression coefficients Pi, P 2 , Pz that are 
given in Eq. (90) by making the normality assumption on the errors of the least squares fit. Then, the 
100(1 — a) percent confidence interval on the coefficients Bj (j = 1,2,3) of the model is defined as (from 
Ref. [4], Eqs. (3.45)): 


CONFIDENCE INTERVALS ON THE REGRESSION COEFFICIENTS 
(... taken from Ref. [4], Eqs. (3.16), (3.17), (3.18), (3.45) . . .) 


Pi - K 


/ 2 , n — p 


a 2 Ci 


— Pj — Pj + t 


:/ 2 , 


a 2 C, 


(91a) 


where Cjj is the jth diagonal element of [X'X] 1 and 
the estimate of the error variance a 2 is defined as follows 


SS Res 

n — p 


1 

n — p 


y' y - P x' y 


(916) 


Equations (91a) and (916) can be used to construct the confidence intervals on the center of gravity 
coordinates after the generic nomenclature used in Ref. [4] is “translated” to the specific nomenclature 
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that is used in the present paper for the calculation of the least squares estimators of the center of gravity 
coordinates. Table 6 below summarizes analogies between the nomenclature of Ref. [4] and the nomenclature 
used in the present paper that were identified after comparing terms of Eq. (32) with terms of Eq. (90). 


Table 6: Comparison of Nomenclature used for Center of Gravity Coordinate Determination. 


Reference [4] 

present paper 

Comments 

Pi, P 2 , p3, 

dx i d y , d z 

regression coefficients; c.g. coordinates 

Pi, P 2 , P3, 

dx ? dy , d z 

least squares estimators; Eqs. (87a), (87b), (87c) 

n 

3 q 

total number of “observations” 

p 

3 

total number of regression model parameter 

*^11) *^12? • • • 5 *^n3 

0 , +7*(1) -w, • ■ • , 0 

coefficient values of model matrix 

Vl, 2 / 2 , 2 In 

Rbal(X), Pbal{ 1), •••, Y ba l{q) 

response variable values 


In addition, because generic vectors and matrices are used in Eq. (916), it is necessary to identify 
corresponding vectors and matrices that were used in the present paper for the calculation of the center of 
gravity coordinates. Table 7 below summarizes analogies between vectors and matrices used in Ref. [4] and 
vectors and matrices used in the present paper that were identified after comparing Eq. (3.13) of Ref. [4] 
with Eqs. (35a), (356), and (376) of the present paper. 


Table 7 : Comparison of Vector and Matrix Nomenclature. 


Reference [4] 

present paper 

Rows X Columns 

Comments 

y 

B 

3 q X 1 

Eq. (35b) 

y' 

b t 

1 x 3? 

transpose matrix 

X 

A 

3 q X 3 

Eq. (35a) 

X' 

A t 

3x3 q 

transpose matrix 

p 

X 

3 x 1 

Eqs. (37b), (87a), (87b), (87c) 

p 

X T 

1 x 3 

transpose matrix 


Now it is possible to identify the remaining terms in Eqs. (91a) and (916). First, Cjj, i.e. , the jth 
diagonal element of the matrix [X'X] -1 needs to be obtained. We know, after looking at Table 7, that the 
following relationship applies: 


X' X = A t A 


Mil U\2 U13 

U2 1 V-22 U23 
U31 U32 U33 


(92 a) 


where the coefficients Un,Ui2, ■ ■ • , 1/33 are defined in Eqs. (82a) to (82/). The unknown coefficients of the 
inverse matrix of [X'X] can be expressed as follows: 


[ x' x p 1 


1 


' C n 

C12 

C13 


= 

C21 

C22 

C23 



1 

C32 

G33 


(926) 


We also know that the product of a square matrix with its inverse equals the identify matrix. Therefore, 
we can write: 


A 1 A 


A 1 A 


= I 


1 0 0 
0 10 
0 0 1 


( 93 ) 
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Each coefficient of the inverse matrix can easily be found by using Cramer’s Rule. It is only required to 
identify the related three 3x3 determinates that have to be constructed by using the coefficients of matrix 
A t A and the coefficients of the identity matrix I. Then, by inspection, we get for Cn, C 22 , and C 33 


C n 

= 

Ui / U 

(94a) 

C 22 

= 

u 2 / u 

(946) 

C 33 

= 

U 3 / u 

(94c) 

U 2 , and C / 3 are defined as 



Mil 

M12 M13 


U = 

M21 

M22 M23 

(95a) 


M 3 1 

M32 M33 



t/i = 


U 2 = 


C /3 = 


1 U 12 

0 U 22 

0 w .32 


u n 
U 21 
U 31 


u 13 

M23 

M33 

U13 

U23 

U 33 


= U22 U33 ~ U23 U-32 


u 11 W12 

M21 U22 

M31 M32 


= Mil M33 


— Mu M22 


M13 M31 


Ml2 M21 


(956) 

(95c) 

(95d) 


and Mn,Mi 2 , ■ ■ • ,M 33 are defined in Eqs. (82a) to (82/). The term y'y — /3 X'y of Eq. (916) still needs to 
be identified before an explicit solution of the confidence intervals on the center of gravity coordinates can 
be assembled. From Table 7 we get: 

y'y = B T B (96) 

We also know after looking at Eqs. (336) and (356) that the following relationship applies: 


B t B = £ T ( 1 ) L(l) + ••• + £ T (g) L(g) 


— -Rbaz(l) + -ffc 2 a ;(l) + Efcai(l) + ••• + R-balil ) + PbaM) + ^bazC*?) (97) 


— [ R-baliQ + PbaliO + ^bai(C) ] 

£=1 

From Table 7 we also know that the following relationship is valid: 

]j X' y = X T A t B 

In addition, after taking the transpose of the left hand side of Eq. (376), we get: 

X [ ci x dy d z ] 

L J 1x3 y 

We also know from Eq. (83) that the following relationship applies: 




Ml 

A t b 

3x1 

M2 



. V 3 . 


(98) 


(99a) 


(996) 
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Finally, after computing the scalar product of the right hand sides of Eqs. (99a) and (996), we get: 

X T A t B = d x ■ v\ + d y ■ V 2 + d z ■ v 3 (100) 

At this point, all generic terms used in Eqs. (91a) and (916) have been identified. Then, after using the 
contents of Table 6 in combination with Eqs. (94a), (946), (94c), (97), and (100), we get the following final 
explicit equations for the confidence intervals on the center of gravity coordinates: 


EXPLICIT SOLUTION OF THE CONFIDENCE INTERVALS 


ON THE CENTER OF GRAVITY COORDINATES 


dx 

t»/2, 39-3 ' G x 

< 

dx — d x + t a / 2 , 3q-3 ' G X 

(101a) 

dy 

^a/2, 3q-3 ’ Gy 

< 

dy 2» dy + t a/2 3q 3 ■ Cy 

(1016) 

dz 

^a/2, 3q-3 ' Gg 

< 

dz d z + t a/2 3 g — 3 • C z 

(101c) 

where the standard errors C x , 

C y , and C z of the coordinates are defined 

as 



c x 

= %/ C u 

(101d) 



Cy 

= \J (T 2 C22 

(lOle) 



c z 

= \J <y 2 C33 

(101/) 


and the estimate of the error variance t? 2 is defined as 



a 2 = 

3 q 

-3-3- ' I G{q) H(q) } 

(1015) 


G{q) = 

Q 

E 

[ Rbal(€) + PbuliO + Y 2 a l{Q ] 

(1016) 



f=i 




H(q) = 

dx ’ 

V\ + dy -v 2 + d z ■ v 3 

(101*) 


Table 8 below lists equations that define the omitted parameters that are used in Eqs. (101a) to (101*). 

Table 8: List of Equation Numbers. 


Parameter 

dx 

dy 

d z 

Cn 

C22 

c 33 

Vl 

V2 

v 3 

Equation 

(87 a) 

(876) 

(87c) 

(94a) 

(946) 

(94 c) 

(84a) 

(846) 

(84c) 


One parameter in Eqs. (101a), (1016), (101c) still needs to be specified before the confidence intervals on 
the center of gravity coordinates can explicitly be calculated. It is the percentage point of the t-distribution 
(i.e., t a/2 39 _ 3 ). This constant depends on (i) the confidence interval choice (e.g., 95 % or 99 % confidence 
interval) and (ii) the number of data points (i.e., q) that are used for the least squares fit associated with 
the center of gravity coordinates. Table 9 below lists the percentage points for different confidence interval 
choices and data point numbers. 
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Table 9: Percentage Points of the t Distribution for p = 3 (from Ref. [4], p. 515, Table A. 3). 


Cl 

a 

a/2 

q 

v = 3q — 3 

^a/2, 3q — 3 

95 % 

0.05 

0.025 

3 

6 

2.447 

95 % 

0.05 

0.025 

5 

12 

2.179 

95 % 

0.05 

0.025 

10 

27 

2.052 

99 % 

0.01 

0.005 

3 

6 

3.707 

99 % 

0.01 

0.005 

5 

12 

3.055 

99 % 

0.01 

0.005 

10 

27 

2.771 


Table 9 shows that at least five data points (i.e., q = 5) are needed in order to obtain percentage points 
of the t-distribution that are near the value of 2 for the 95 % confidence interval and near the value of 3 
for the 99 % confidence interval. Therefore, the authors recommend to use at least five data points for the 
least squares fit that is associated with the regression analysis associated with the determination of both the 
model weight and center of gravity coordinates. 
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Appendix 5: General Error Analysis of the Model Weight Prediction 

An estimate of the upper bound of the wind tunnel model weight prediction error was derived in App. 2 
of the paper. A more general error analysis of the model weight prediction can be performed by applying a 
well known error propagation formula from the literature to Eq. (30). This formula estimates the variance 
in the model weight prediction (a discussion of the error propagation formula is given, e.g., in Ref. [8], p. 49). 

The general analysis of the model weight prediction error and the application of the error propagation 
formula starts by expressing Eq. (30) in a slightly different format. The least squares estimator of the model 
weight can be described by the following set of equations: 

i « 

w = — • u/(£) (102a) 

q l=i 

where 

™'(0 = y/AUS) + SUO + N 2 al (0 (102 b) 

It can be concluded from the above equations that the final model weight estimate is a function of the 
individual axial force, side force, and normal force measurements of each data point that is used to compute 
the model weight estimate. This conclusion can be summarized as follows: 

ui = F[A ba i{ 1), ... , A ba i(q), S bal (l), ... , S ba i(q ), N ba i( 1), ... , N bal (q) ] (103) 

Now, it is possible to apply the error propagation formula from the literature to Eq. (103). The error 
propagation formula is defined for a multivariate function Y where X \ , . . . , X $ are a set of independent 
variables. This characteristic of function Y can be summarized as follows: 

Y = F{ X u X 2 , ... , X+] (104a) 

The error propagation formula is given by the following formula (taken from Ref. [8], p. 49, Eq. (3.3)): 

^ y > = (w) 2 '" 2 ™ + - + (w) 2 '" 2 ™ (1044) 

The error propagation formula can be applied to Eq. (103) after comparing the nomenclature used in 
Eq. (103) with the nomenclature used in Eq. (104a). Table 10 below lists corresponding symbols. 

Table 10: Comparison of Nomenclature used by the Error Propagation Formula. 


Nomenclature used in Eq. (104a) 

Nomenclature used in Eq. (103) 

Y 

W 

X±, X 2 , , X+ 

Af,ai(l)> • • 

• j A-baliSl) i ••• ? ^bal^St) i •• 

, N ba i(q) 

$ 

3 q 


Then, after replacing the generic symbols used Eq. (104&) with the symbols listed in the second column 
of Table 10, we get the following error propagation formula for the model weight estimate 


< ’ 2( ” > = ■ < ’ 2( ' 4 “ (1)) + - + (a^rw) • " 2(Ar “ (,)) (105) 

where a 2 (w ), a 2 (A ba i), cr 2 (S ba i), and cr 2 (N ba i) are the variances in the model weight, axial force, side force, 
and normal force predictions. It is reasonable to assume that the variances in the axial force prediction are 
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identical for all “wind-off” data points. The same assumption can be applied to the side force and normal 
force predictions. These conclusions can be expressed as follows: 


a(A ba i(l)) 

= cr(A ba i(2)) = ••• 

= cr(A ba i(q)) = 

const. = 

a(A) 

(106a) 

&(S ba i{ 1)) 

= v(Sbai{2)) = ••• 

= <r(S ba i(q)) = 

const. = 

cr(S) 

(1066) 

<r(iWl)) 

= a{N bal { 2)) = ••• 

= cr(Nbai(q)) = 

const. = 

a(N) 

(106c) 


Partial derivatives used in Eq. (105) are obtained after applying the chain rule to Eq. (102a) and (1026). 
Then, we get: 


d w _ A ba i(£) 

d A ba i(£) q-W(C) 

d W _ S b al(£) 

d S ba i((i) q-w'if) 

d W _ N ba i (£) 

d N ba i(£) q-w'{£) 


(107a) 

(1076) 

(107c) 


Finally, after using the right hand sides of Eqs. (106a), . . . , (107c) to replace related terms in Eq. (105), 
we get the following relationship for the variance in the model weight prediction: 


er 2 (w;) 


VARIANCE IN THE MODEL WEIGHT PREDICTION 



(108) 


The formula of the variance in the model weight prediction, i.e., Eq. (108), can significantly be simplified 
by taking subject matter knowledge about both wind tunnel test procedures and strain-gage balances into 
account. The roll angle of the wind tunnel model, for example, is often either 0° or 180° when “wind-off” 
data points are recorded. Consequently, the gravitational acceleration has no component in the direction of 
the side force. This conclusion can be summarized as follows: 

S b ai(0 « 0 (109a) 

We also know that the contribution of the normal force to the model weight estimate is significantly 
larger than the contribution of the axial force because (i) the primary sensitivity of the normal force of a 
typical strain-gage balance is one order of magnitude larger that the primary sensitivity of the axial force 
and (ii) the pitch angle of the wind tunnel model is in most cases less than 20° when “wind-off” data points 
are recorded (see also the related discussion in App. 2). Therefore, the following inequality applies: 

I N bal \ > \A bal \ for |/i| < 20° (1096) 

Consequently, the weight estimate u/(£) defined in Eq. (1026) for each individual “wind-off” data point 
with index £ can be approximated by using the relationship: 

«>'(0 = y/AUO + SUZ) + NLtt) =► ( ,A =► ^(O « (109c) 

[ A bal \ sj ^ A bal \ sJ J 

In addition, the variance in the normal force measurement is expected to be significantly larger than the 
variance in the axial force measurement because, as mentioned earlier, the primary sensitivity of the normal 
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force of a typical strain-gage balance is one order of magnitude larger that the primary sensitivity of the 
axial force. This conclusion can be summarized as follows: 

cr(iV) er(A) (109d) 

Then, after using the relationships given in Eqs. (109a), . . . , (109d) in Eq. (108), the following approx- 
imation of the variance in the model weight prediction is obtained: 

a 2 {w) « • ^V(TV) (110) 

q l=i 

Equation (110) can be simplified further because we know from Eq. (106c) that the variance in the 
normal force measurement is a constant. Therefore, the following relationship is valid: 

J2<7 2 (N) = q ■ ct 2 (N) (111) 

5=1 

Finally, after using the right hand side of Eq. (Ill) to replace the summation term on the right hand 
side of Eq. (110), we get the following approximation of the variance in the model weight prediction: 


APPROXIMATION OF THE VARIANCE IN THE MODEL WEIGHT PREDICTION 


a 2 (w ) 


° 2 {N) 

q 


( 112 ) 


Equation (112) can be used to estimate the minimum number of “wind-off” data points needed to satisfy 
a model weight prediction tolerance TOL(w ) that a wind tunnel customer may specify. The corresponding 
relationship can be obtained by first taking the square root on both sides of Eq. (112). Then, we get: 


cr(w) 


°(N) 

VT 


(113) 


By inspection, the following inequality has to be satisfied in order to meet a customer’s tolerance for 
the model weight prediction: 

- - < TOL{w) (114) 

Finally, after solving inequality (114) for the total number of “wind-off” data points, we get the relationship: 


MINIMUM NUMBER OF “WIND OFF” DATA POINTS 


° 2 {N) 

TOL 2 (w) 


(115) 


Inequality (115) clearly shows that the minimum number of “wind-off” data points needed to meet a 
customer’s tolerance is directly proportional to the variance in the normal force measurement. This variance 
in the normal force measurement is not simply the square of the standard deviation of the back-computed 
normal force residuals that are obtained after completion of a balance calibration data analysis. It also 
includes errors associated with test procedures, instrumentation, and hardware that are used to obtain the 
“wind-off” data points in the test section of the wind tunnel. 


33 

American Institute of Aeronautics and Astronautics 


Gage 

Output 



A 


output change used 
for load calculation 




global output datum = natural zero 


o 



Pitch Angle 


Fig. 1 Balance gage output parts of a “wind-on” data point. 



z 



Fig. 2 Definition of Cartesian coordinates and load signs in the “primary” balance axis system. 
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z 



Fig. 3 Description of the gravitational acceleration and the center of gravity co- 
ordinates of a wind tunnel model in the “primary” balance axis system. 
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Gage 

Output 


“ wind-off ” 



output change used 
for load calculation 


Pitch Angle 


Fig. 4 Balance gage output parts of a “wind-off” data point. 



Fig. 5 Definition of Cartesian coordinates and load signs in the “alternate” balance axis system 
(the definition of the coordinates and load signs is copied from Ref. [1], p. 3, Fig. 1). 
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“Wind— On” Test Condition 

“Wind Off” Test Condition 

Comments 

Balance Loads of a 

Balance Loads of Point Set 

Loads are obtained after the appli 

Single Data Point 


cation of the load iteration process 


Abal(£), Sbal(£)i ■■■ , YbalU) 

that uses (i) the difference between 

-A-bal 5 ^ bal > ■ • ■ j ^ bal 


raw outputs and natural zeros and 


where 

(ii) the data reduction matrix as 
input; “wind-on” and “wind-off” 


VI 

VI 

i-H 

loads are computed using the same 
natural zero set for the iterations. 

Aerodynamic Loads 

Aerodynamic Loads 

The aerodynamic loads on the model 

of a Single Data Point 

of Point Set 

in the “wind-on” test condition 
are a primary result of the wind 

A-aero — bal A weight 

Aaero (0 = 0 

tunnel test. 

^aero — ^ bal & weight 

Saero( 0 = 0 


Y a ero — ^ bal ^ weight 

Yaeroit) = 0 


Unit Gravity Vector 

Unit Gravity Vector 

Vector components are the result 

Components of a 
Single Data Point 

lx, lyi Iz 

Components of Point Set 

ix{£), iy(Oi iz {0 ; i < £ < q 

of model attitude measurements. 

Forces caused by 

Least Squares Estimators 

The least squares estimators of 

the Model Weight 

of the Model Weight and the 

the model weight and the center 


Center of Gravity Coordinates 

of gravity coordinates, i.e., w, d x , 

A weight V x ’ ^ 


d y , d z , and the unit gravity vector 


w ~ w => Eq. (30) 

components of the “wind-on” test 

S weight = Vy ’ ^ 


condition, i.e., 7 X , j y , 7 z , are used 


d x « d x =>• Eq. (87a) 

as input for the calculation of 

N weight = ^fz * ^ 


forces and moments associated with 


d y ss d y => Eq. (876) 

the model weight; the least squares 

Moments caused by 


estimators are constant as long as 

the Model Weight 

d z S3 d z =>■ Eq. (87c) 

the weight distribution and the 
total model weight do not change. 

Rweight = lz- W ■ dy- "1y ■ W ■ d z 

Confidence Intervals on 
the Model Weight and the 


Pweight = Tx ’ ^ ’ d z ^fz ’ W ' d x 

Center of Gravity Coordinates 


Y weight = Vy ' ^ * d x ^fx ‘ W ’ dy 

w =£* Eqs. (55a) to (55d) 
d x , d y , d z =>• Eqs. (101a) to (101*) 



Fig. 6 Summary of basic elements of the new method to predict model weight corrections. 
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